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We show that the kinetics of microribbon formation of amelogenin molecules is well described 
by a combination of translational and rotational diffusion of a simplified anisotropic bipolar model 
consisting of hydrophobic spherical colloid particles and a point charge located on each particle 
surface. The colloid particles interact via a standard depletion attraction while the point charges 
interact through a screened Coulomb repulsion. We study the kinetics via a Brownian dynamics 
simulation of both translational and rotational motions and show that the anisotropy brought in 
by the charge dramatically affects the kinetic pathway of cluster formation and our simple model 
captures the main features of the experimental observations. 
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The self-assembly of colloidal particles and globular 
proteins with inhomogenous properties into various de- 
sired structures is an important area of current research, 
in which the inhomogeneity (such as an amphiphilic 
group and a residual group) plays a crucial role in the 
kinetic pathway of cluster formation and eventually im- 
pacts on the stabilized crystal structures. Examples in- 
clude Janus particles and DNA-coated colloidal particles 
that lead to the next generation of building blocks of new 
materials and have potential applications in fabricating 
phontonic crystals, targeted drug delivery, and electronic 
equipment Another important example is amelo- 

genin, the chief hydrophobic protein in enamel matrix 
with a hydrophilic 2 5- amino acid C terminus. Amel- 
ogenin is involved in the mineral deposition and is re- 
sponsible for the major structural process during enamel 
biomineralization. A recent experiment shows that the 
hydrophilic C-terminus is essential for the self-assembly 
of amelogenin into microribbons, which may also be rel- 
evant to the elongated and oriented growth of apatite 
crystals during biomineralization 0] . The experiment 
also revealed that the microribbon is formed by a hierar- 
chical organization of amelogenin molecules into a chain 
of nanospheres. 

In this article we present a coarse-grained model 
of amelogenin and explain how the microribbon forms 
through self-assembly. The amelogenin molecule is hy- 
drophobic with a charged hydrophilic tail. Our model de- 
scribes this in a simplified way, representing the monomer 
as a spherical molecule, with the charged hydrophilic tail 
replaced by a single tethered point charge located on the 
surface of the molecule. Using Brownian dynamics simu- 
lations, we investigate the static and dynamic properties 
of the self-assembly process. We show that the anisotropy 
brought by the charge dramatically affects the kinetic 
pathway of cluster formation and our simple model cap- 



tures the main features of the experimental observations 
in the early stages. 

In the experimental studies of amelogenin assembly 
process, one typically adds salt and a precipitant such 
as polyethylene glycol (PEG). Although PEG-protein in- 
teractions are quite complicated 0, we model this as 
a depletion interaction that we take to be given by the 
Asakura-Oosawa potential plus a repulsive hard-core-like 
interaction, depending on the center-center distance be- 
tween spherical particles, Udri^) = UAo{rij) + Uhdr'^j), 
where 
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55 The cut-off range Tc = 1 + ^, where ^ is the size-ratio 

56 between a PEG coil and a colloidal amelogenin particle 

57 that controls the range of the depletion interaction, and 

58 4>p is the value of the PEG volume fraction that controls 

59 the strength of the interaction described by the absolute 

60 value of the minimum potential depth = \Umin\ 

61 The hard core potential is given by 



kT 
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62 We set a = 36, since the values of a < 36 have been 

63 reported to lead to anomalies when a mimic of the hard- 

64 core potential is required in the potential [IQ] . The 

65 point charges interact with each other through a screened 

66 Coulomb potential, 
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in which magnitude controlled by e and range controlled 
by Debye screening length Ad, are varied in the simu- 
lations. We find that the early morphology of the self- 
assembly is more sensitive to the size of A^i rather than 
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71 e. These Coulomb charges exert torques on adjacent 

72 molecules and hence produce a rotational motion of the 

73 molecules that is included in our Brownian equations of 

74 motion, read as: 



mf, = -V(C/^ + Uf) - Tth + W,{t), 



luji ^ Ti - TrUJi 



(4) 
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76 where m, /, 7\, uji, Ti are the mass, moment of iner- 

77 tia, position vector, angular velocity, and torque, respec- 

78 tively, of the ith colloidal particle. The mass of the point 

79 charge is ignored in this model. (r^) is the transla- 

80 tional (rotational) friction coefficient. Wi and are 

81 the random forces and torques acting on the ith col- 

82 loidal particle respectively, which satisfy a fluctuation- 

83 dissipation relation {Wi{t) ■ Wj{t')) = &kTTt5ij5{t - t'), 

84 • M>'(t')) = QkTTr5rj5{t - t') We choose 

85 Tt = 0.5, Tr = 0.167 and the time step At = 0.005 in 

86 reduced time units of a{m/kTYl'^ with m = 1 [l3|. The 

87 other parameters are chosen according to the range of 

88 values of the experimental data displayed in Du et al's 

89 paper Q. In particular, here the results shown are for 

90 the values C=0-1, Um = 6kT, Xd — 0.4 (ah length scales 

91 are measured in units of monomer diameter and energies 

92 are scaled by e). We consider a small volume fraction of 

93 the molecules, typically of the order of / = 0.01 or 0.02. 

94 Periodic boundary conditions are enforced to minimize 

95 wall effects. All simulations start from a random initial 

96 monomer conformation and the results for the kinetics 

97 are averaged over several (5-10) runs. 

98 We show typical configurations in the early stage of 

99 microribbon formation in Fig. 1. We see that ini- 

100 tially the monomers form oligomers and the oligomers 

101 self-assemble into larger aggregates contained (approxi- 

102 mately) within nanospheres, then nanospheres associate 

103 together into a chain structure, quite similar to the pro- 

104 cess of amelogenin self-assembly (see Fig. 3 in reference 

106 i). 

107 The Coulomb forces cause the majority of the charges 

108 to point outwards inside the nanospheres, as shown in 

109 Fig. 2. This again is reminiscent of the situation 
no with the hydrophilic tails of amelogenin. To character- 
in ize the orientational distribution of point charges within 
112 a nanosphere, we introduce a parameter 6, defined as 

I y, I I , where u„ is the position vector of 

a point charge referenced to the center of its host colloid 

115 and Uc the position vector of a colloid referenced to the 

116 center of mass of this nanosphere. Therefore 6 < 90° 

117 if the charge points outwards and 9 > 90° if it points 

118 inwards. Due to fluctuations and our choice of short re- 

119 pulsion interaction range, the distribution of 9 is spread 

120 around a small nonzero peak position, as shown in Fig. 

121 2, for the configuration in inset. Subsequent aggregation 

122 forms via necks between nanospheres, leading to larger 

123 structures that eventually form (flexible) microribbons. 
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FIG. 1. (color online), (a) The model of amelogenin molecule 
consists of a spherical colloid particle and a point charge lo- 
cated on its surface that preserves a bipolar nature, (b) and 
(c) Oligomerization of amelogenin molecules occurs by means 
of hydrophobic interactions and modification due to Coulomb 
repulsions, (d) Nanosphere structures are formed through 
aggregation of monomers and oligomers, (e) Further associa- 
tion of nanospheres results in larger assemblies among which 
chains are formed. 




FIG. 2. (color online). Distribution of 9; 9 shows whether the 
point charge orients outwards (9 < 90°) or inwards (6^ > 90°). 
Inset: Nanosphere structure, large (grey) spheres indicate 
amelogenin molecules while small (red) ones for point charges. 



125 We characterize the morphology of the clusters in 

126 terms of their fractal dimension Df. The q-dependence 

127 of S{q,t) shown in Fig. 3(a) is given by a slope of - 

128 2 on a log-log plot over a reasonable range of q. This 

129 value oi Df is larger than the typical diffusion-limited- 
no cluster-cluster aggregation (DLCA) value of 1.8 as the 

131 ribbon-like clusters observed here have significant short 

132 range ordering [l^. For this reason we speculate that 

133 the repulsive Coulomb force causes the surface particles 

134 to reorganize, which reduces the surface roughness of our 

135 model, in contrast to the case without the Coulomb in- 

136 teractions (see reference [l3l)- We find that a peak de- 

137 velops in the structure factor whose magnitude increases 

138 with time. Correspondingly, the peak position decreases 

139 as a function of wave number with increasing time; these 

140 features are typical of spinodal decomposition 14]. How- 
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141 ever, we note that for deep quenches that lead to cluster 

142 growth considered in this model, the system is controlled 

143 by two characteristic lengths that evolve differently in 

144 time and an apparent scaling for the structure factor can 

145 only be observed over some period of time when these 

146 two characteristic length scales become comparable to 

147 each other 
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FIG. 3. (color online). Log- log plot of structure factor at 
several times. The dashed line indicates fractal clusters with 
S{q) (X q~^f , where Df — 2.0, while the dotted line indicates 
the Porod regime S{q) oc g"'''"'"^', d — 3. 



150 The kinetics of the cluster growth process in this sim- 

151 pie model for amelogenin self-assembly is consistent with 

152 a cluster-cluster aggregation mechanism, as shown in Fig. 

153 4. The number of clusters decreases inversely with time 

154 (i.e. the kinetic exponent ^ = 1) and the radius of gyra- 

155 tion increases as a power law with an exponent of n = 0.5 

156 17|, |18| . The relation between z and n involves the frac- 

157 tal dimension in the following way: n — z/Df. Thus the 

158 kinetic exponents are consistent with a fractal dimension 

159 of I?/ = 2 as well. 




FIG. 4. (color online). Plot of (a) number of clusters, Nc and 
(b) radius of gyration, Rg as a function of time in a log-log 
plot. The dotted line in (a) has a slope of —1; the dashed line 
in (b) has a slope of 0.50. 



In order to check the stability of the clusters produced 
in straight simulation from a quench (Um = 6kT) into the 
two phase gas-solid region, we "heat" the system formed 
after 5000 steps from the initial quench Um = QkT to 
Um — '^kT and run the simulation for another 5000 time 
steps. This heating corresponds to the change in PEG 
concentration as it is done in the experiment. A typi- 
cal run is shown in Fig. 5. As shown in the figure, we 
find that the microribbons arc quite robust under such a 
change of conditions, as found experimentally for a wide 
range of conditions (including a change in the PEG con- 
centration). We have also checked the structure factor 
behavior for the heating process and found that the struc- 
ture factor curve has the same exponent of -2 as at the 
cooler temperature (see Fig. 3). This is consistent with 
the robustness of the major cluster morphology. 





FIG. 5. (color online), (a) Morphology of entire system (^ = 
0.1, / = 0.02) from another 5000 time steps shallow quench 
of Um ~ 4fcr after first quench into Um ~ QkT, 5000 time 
steps, (periodic boundary condition are enforced), (b) One 
of the clusters inside (a). 



180 We characterize these microribbons by computing 

181 three eigenvalues of the gyration tensor, denoted as , 

182 A|, and A| in descending order [l^l- We calculate the ra- 

183 tios among these three principal values, which character- 

184 ize the relative colloid (or charge) particle distribution. A 

185 typical result is given in Table 1, for the colloid particles 
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TABLE I. Typical ratios of eigenvalues of gyration tensor 







A2/A3 


Colloid particles 


9.57 


9.03 


Point charges 


9.25 


8.74 



186 and point charges separately, based on the configuration 

187 in Fig. 5(b). These ratios of eigenvalues demonstrate 

188 the anisotropy of these clusters, suggesting that the clus- 

189 ter is ribbon-like. The difference of the ratios between 

190 colloidal particles and point charges may indicate that 

191 under electrical repulsion, the rotation of monomers in- 

192 side the cluster leads to the small reduction of anisotropy 
im in the point charges' spatial distribution. 

195 From the gyration tensor we can determine the typical 

196 length scales of these clusters by taking square roots of 

197 the eigenvalues. To illustrate this, we use the colloid con- 

198 figuration shown in Fig. 5(b), which gives 18.27, 17.75, 

199 and 5.91 for the three length scales. Since the chain struc- 

200 ture is the association of several nanospheres (as shown 

201 in Fig. 1), the typical diameter of the nanosphere in 

202 this configuration is at the most 5.91. We can also use 

203 this estimate of the diameter to obtain an estimate of 

204 the minimum wavenumber that corresponds in Fig. 3 

205 to the crossover to the Porod regime (corresponding to 

206 compact clusters on a local scale). This yields a value 

207 of Qmin — 27r/5.91 = 1.06, which is consistent with the 

208 behavior shown in Fig. 3(b). Since the hydrodynamic 

209 radius Rh of the amelogenin macromolecule is about 2.2 

210 nm, the Rh of the nanosphere for the configuration in 

211 Fig. 5(b) is 2.2 x 5.91 = 13.00 nm; this lies in the range 

212 of values of the radii of typical nanospheres, 10 ^ 25 nm, 

213 found in the experiment (cf. Fig. 3, reference Q). 

214 In summary, we have developed an anisotropic, bipo- 

215 lar model for the hierarchical self-assembly of amelogenin 

216 molecules and have carried out Brownian dynamics sim- 

217 ulations of the self-assembly process. Simulations show 

218 a hierarchical self-assembly process where the molecules 

219 aggregate to form dimers, hexamers to nanospheres, and 

220 the assembly of the nanospheres then lead to the for- 

221 mation of microribbons in agreement with experimental 

222 findings. The relative strengths of the interaction param- 

223 eters can lead to a phase diagram where the nanospheres 



224 are stable against microribbon formation. Thus a con- 

225 trol over the morphology of the clusters can be achieved. 

226 Such a study is currently underway. 
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